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We consider a process of the Andreev reflection between a normal metal and the s-wave super- 
conductor in the FFLO state. It is shown that the process takes place if the energy of the incoming 
electron is bound within the finite interval called the Andreev window. The position of the window 
determines the value of the non-zero total momentum of Cooper pairs and the value of the gap. 



During the last two decades one has faced a renaissance of the interest in non-uniform superconductivity, where 
the spatial symmetry is broken by a non-zero total momentum q of Cooper pairs [l[ . This has happened mainly due 
to the experimental discoveries of the possible candidates for the superconductors in the FFLO (Fulde-Ferrel-Larkin- 
Ovchinnikov) state Q. Recently, strong evidence for the existence of the FFLO state has been brought forward for 
the organic superconductor Q and for the CeCoIn5 superconductor Additionally a new interest has arisen in 
the field of strong interaction physics where the new state of matter called color superconductivity was suggested [f| • 
In condensed matter physics the FFLO state requires an applied external magnetic field which leads to the Zeeman 
splitting of the Fermi surface of conduction electrons, while in QCD the Fermi surfaces of quarks are already split 
because of the different masses of different flavor quarks @ . It is also a good place to point out that the non- uniform 
condensates were also considered in the context of chiral symmetry breaking [7|, the subject which is still under the 
current debate Q. 

The Andreev reflection 0] between a normal metal and an anisotropic superconductor with a directionally dependent 
gap was discussed in (lp| with the main interest concentrated on the d-wave superconductors. In this paper we consider 
the junction between a normal metal and the s-wave superconductor in the simplest possible FFLO state where the 
gap parameter has an oscillating phase A(r) ~ exp (iq ■ r). All the calculations are performed in the pure Pauli limit. 
We show that the Andreev reflection can be used, at least in principle, to detect the existence of the non-zero total 
momentum of Cooper pairs. Despite the clear result still a lot of work must be done to bridge between the ideal case 
we discuss in this paper and the description of real materials used in experiments. 

In the first section we discuss the Bogolubov - de Gennes equations and the general properties of their solutions. 
In the second section an exact, numerical solution is given for some generic parameters. In the last section there is a 
short discussion of the possible further work. 



Bogolubov - de Gennes equation for a non-relativistic superconductor in one dimension takes the form: 



where the FFLO gap function is chosen as A(z) = A e 2l<?z and the FFLO momentum q is assumed to be non-negative. 
The plane wave ansatz f(z) = fk+ q e l ^ k+q ^ z and g(z) — gk-q^ k ^ q ^ z lead to the matrix equation 



I. INTRODUCTION 



II. BOGOLUBOV - DE GENNES EQUATIONS FOR THE 1-DIM LOFF PHASE 
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which gives the dispersion relations 

E ± = \ (tk+ q - £k- q ± v /(e fc _ 9 + e fe+9 ) 2 + 4|A | 2 ^ . (4) 

In Fig 1, examples of the dispersion relation ((4]) are given for q = (left panel) and q ^ (right panel). It can 
be immediately seen that the non-zero value of the total momentum q breaks the parity symmetry k — > —k. The 
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exact shape of the curves depends on the size of the total momentum q. For the values of < E q < 1 4 ° 1 , where 

2 

E q = there are two positive minima of E + of different depths at the points that are not related by the parity 
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transformation. The curve E- resides below zero. However, for E q > 1 4 ° 1 one of the minima of E+ descents below 
zero. Simultaneously a part of the curve EL emerges above zero on the opposite side of the momentum axis. This 
phenomenon leads to more complicated dispersion relations for hole-like and particle- like excitations. 

Let us consider the process of the Andreev reflection in a simple one dimensional geometry configuration where the 
FFLO superconductor resides at z > with the gap parameter A(z) = Q(z)Aoe 2lqz . The normal metal is at z < 
with a junction to superconductor at z — 0. An electron of the mass m and the energy E above the Fermi energy 
arrives from the left at the junction. 




FIG. 1: The dispersion relations Q with the upper curve E+ and the lower curve E-. In the left panel q — and in the right 
panel q 7^ 0. The lines a, b represents the energy E of the incoming electron. The crossing points k^ 2 describes the solutions 
of the equations E = E± . 



Then one needs to solve the equations ([1]) for fk+ q and gk-q with the sewing conditions ^><(0) = ^>>(0), ^<(0) = 
ip>(0) imposed at z — 0. The incoming electron of given energy E can excite quasiparticles inside the superconductor 
with real or complex momenta. Quasiparticles with the real momenta can propagate inside the superconductor freely 
whereas those with the complex momenta penetrate only the region close to the junction. The solutions with real 
momenta exist if there are non-zero real solutions of the equations E± = E. Let us first consider the case of normal 
superconductor with q = 0. These solutions geometrically are placed at the crossing of the horizontal lines a, b and 
the curves E± (lines a,b in the left panel of Fig. 1). For the s-wave BCS superconductor there are four possible real 
solutions for E > A (line b) among which kf describes hole-like and k 2 particle-like excitations propagating from 
the left to the right whereas fcf , k 2 are those which propagate in the opposite direction. Only two of these solutions 
are compatible with the boundary conditions at infinity (determined by the direction of the incoming electron plane 
wave). Let us consider the case where the electron is incoming from the left. Then solutions fcj~ 2 are rejected and the 
only left are kf 2 - It is worth to mention that the incoming electron mainly excites the electron- like quasiparticle in 
the superconductor, which follows from the kinematics of the reflection process. In the case of E < A there are no 
real solutions (no crossing of the line a with the curve E + in Fig. 1, left panel). The lack of the freely propagating 
excitations in the superconductor results in the process called the Andreev reflection. 

In the case of the FFLO superconductor one can immediately infer from Fig. 1 (right panel) that there are always 
real solutions of the equations E± = E. Let us consider first the solutions described by the horizontal line a. This 
corresponds to the case of E > A for the uniform BCS superconductor. There are two solutions k^ 2 of possible 
four that are compliant with the boundary conditions. Similarly to the uniform case the electron-like excitation k 2 
dominates over the hole-like k^ quasiparticle. In the case described by the line b there are two possible solutions 
fej , k 2 and only one solution k^ is compatible with the boundary condition. This is exactly the one which describes 
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the hole-like quasiparticle propagating from the left to the right. However, this solution is strongly suppressed by 
the kinematics of the reflection. As a result the quasi- Andreev process takes place when the energy of the incoming 
electron is within the gap between the minimum of the upper E + and the maximum of the lower E- . Approximately 
this is in a range 

^-^2m^ -q 2 - |A | < E < —y/2mn - q 2 + |A | (5) 
m m 

with the corrections of the order of ^° ^ ne um it s - These inequalities are valid as long as E q < Ep. 

The inequalities ([5]) in principle can be used for the determination of the important superconductor parameters. 
The difference between the upper and the lower bound gives the gap parameter |A | whereas the sum determines the 
value of the total pair momentum q when the Fermi energy /i and the mass of the charge carrier are given. 

III. NUMERICAL RESULTS 

We consider an electron that is injected from the conductor side of the junction. In this case the wave function 
takes the form 

M*)= (l) eikZ + B ( jV^ + C^e^, z<0, (6) 



*>«-(J:)-K^)^' + '(^)^' >a (7) 

The total probability current obeys a continuity equation 

ip = i (Vl^i +n^h -^:) ■ (8) 

One can decompose the current j p on both sides of the junction into parts connected with (quasi)hole and 
(quasi)particle excitations. These are the incident probability current generated by the incoming electron jf = ^-k, 
the probability current connected with the reflected hole jf h = —^P |C| 2 , the probability current connected with the 
reflected electron j< e = — — k \B\ 2 . At the superconducting side one defines the probability current connected with 
the transmitted quasiparticle 

3l = -k| 2 ((l/ 2 | 2 -| 5 2| 2 )5R[fc 2 + ]+9), (9) 
and the probability current connected with the transmitted quasiholc 

j& = ^^((l/il'-lflil'Wl+s)- (10) 

Then the hole reflection coefficient is defined as Rhoie = UwJ/bT"! anc ^ ^ ne quasiparticle transmission coefficient is 
defined as T quasi = |j t >|/|j<|. 



For the presentation of our numerical results we set the typical values of the superconducting gap |Ao| = 0.001 eV 
and the Fermi energy fi — 1 eV. The value of the momentum q corresponds to the wavelength of several interatomic 
distances ~ 10 nm which gives q ~ 10 eV. The quasi- Andereev reflection, T quas i «0 and Rhoie ~lj is expected to occur 
within the energy range (O, that for established parameters induces 0.0187826 eV < E < 0.0207826 eV. As can be seen 
in Fig. 2, this expectation is very well confirmed by the numerical findings. For the incoming electron of the energy 
outside the Andreev window ([5]) there are solutions of the propagating quasiparticles in the FFLO superconductor, 
and they are responsible for the non-zero transmission coefficient T quas i. It is worth to mention that in the opposition 
to the situation with q — 0, there are always non-evanescenting waves in the FFLO. Even within the Andreev window 
([5]) the solutions with the real momentum k exists. The amplitude F (responsible for the propagation of quasiholcs 
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FIG. 2: The probability of the hole reflection (left panel) and the probability of the quasiparticle transition (right panel) as a 
function of the energy (eV) of the incident particle. 



in the FFLO) is non zero in this region and the total transmission coefficient T total = \jte + 3th\l\3i~\ i s different from 
zero. However, the value of \F\ 2 is negligibly small and is not visible in Fig. 4. 

Another interesting quantity is a c harg e transport on the superconductor side of the junction. The total charge 
current obeys the continuity equation [ill ] 



d 



d 



e g- t (\^\ ~ IVU I 2 ) + J? = 4e3[A^%] 



(11) 



where 



<:)z 



-(s^v^H^tW-i]), 

-4e3[A^|]. 



(12) 
(13) 



The quasiparticle charge current j> and the charge current carried by the condensate jf c take the form: 



2 - \n,\ 2 \) P -^( k t> 



j> = -(\Fr(^kt}+q(\hr-\gxn)+\jr(mt}+Q(\M-\9 2 n) 

+ Wf2F*f* 1 (k+* + k+ + 2q)e i(k 2- k t'>-] + $t[Jg 2 F*g* 1 {k+* + fc+ - 2q) e ^ + ^>]), 



4eA S| 



23(fc+) li ' (kt-kt*) 



(i _ e i ( k t- k V) z ^ j 



(k 2 ^ ) 



(14) 



(15) 



The main contribution to the charge current originates from the part of the equation (|15p that is proportional to the 
coefficient I J| 2 , that is: 



2eA |J| 3(/ 2 *g 2 ) (, „-2S(fe+)* 



1 - e 



(16) 



Considering the limit q — > the charge current has a simple form 



j>(q = 0) = 2ev F 1 - e 



1\J\ A| 2 -E 2 



(17) 



where vf — \ — the Fermi velocity. For g ^ one can rewrite equation (|16[) in the approximate form: 



2J|A|2-(^-^V2m/x-g 2 ) 2 
1 — exp . z 



■ r, 1 r 



(18) 



nffi- 0.005 



0.0002 0.0004 0.0006 0.0008 0.001 eV 



0.019 0.0195 0.0200 0.0205 eV 



FIG. 3: The plot of l/£ [nm ] as a function of the energy [eV] of the incoming particle. Left chart q=0, right chart q=10 eV. 



Here one introduces a standard parameter £ that describes the penetration depth of the evanescenting currents 

e(g=o) = - h Z = , e(?) = 



2y/\A\*-E*- 



(19) 



2 v /|A|2-(^-X v /2m/i-^)2 
In Fig. 3 we present the penetration depth £ as a function of energy of the incoming electron. We keep the values 
Jsc,c/2ei' F 1.00 




0.75 - 



50 



0,25 



z [nm] 



FIG. 4: The charge current and the current carried by the condensate in a units of 2evF as a function of a distance z [nm] from 
the junction at z = (black points refers to the incident energy E = 0.02 eV, gray points to energy E — 0.0189 eV ). 

of m, A and /i constant. As one could expect, when q ^ penetration distance approaches towards infinity as E 
approaches towards the lower or higher limit at the Andreev window The minimum value is at E = ^--y/2rn^ — q 2 . 
On the contrary, for q — 0, the maximum penetration distance is at E — A and the minimum is reached at zero 
energy. The dependence of £ on the energy is also seen in Fig. 4, that presents the conversion of the normal charge 
current into the supercurrent at the FFLO phase. 



IV. CONCLUSIONS 



In this paper we have considered a process of the Andreev reflection between a normal metal and s-wave supercon- 
ductor in the FFLO state with a single plane wave in the pure Pauli limit. We have found that the Andreev process 
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takes place only for electrons with the energy located within the Andreev window given by equation This is 
clearly visible in the dependence of transmission and reflection coefficients associated with the energy of the incoming 
electrons. The other interesting quantities are the charge currents flowing through the junction and the penetration 
depth of the excitation inside superconductor. This last parameter depends on the energy of the incoming electron 
in very different way compared to the s-wave BCS superconductor (Fig. 3). In conclusion the Andreev process, in 
principle, can serve as a good probe to look for the FFLO state. Obviously the model given here is a far reaching 
idealization and additional works must be done to bring the results closer to experimental conditions. However, the 
results presented in this paper are robust in a sense that they are the simple consequences of the existence of the 
non-zero total momentum of Cooper pairs. 

The next interesting step is analysis of the FFLO state with an oscillating value of the gap parameter A ~ cos(g-x). 
Another important point is to adopt the external magnetic field in the picture and try to go off the pure Pauli limit. 
One can also extend the analysis of the Andreev reflection in the color superconductors [12| for the FFLO states. 
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